Observation of the (j) —> 7/7 decay with four 
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Abstract 



Using 11.6 pb of data collected in the energy range 0.864-1.06 
GeV by CMD-2 at VEPP-2M, the rare decay mode 4> -> rj'j was 
observed via the decay chain rf — > it + it~7], r\ — > 7r + 7r _ 7r° or rj — > 
7r + 7r _ 7. The following branching ratio was obtained: 

Br(<f> -> 7/7) = (49li j ± 0.6) • 10~ 5 . 

1 Introduction 

Measurements of radiative decays of light vector mesons provide good tests of 
the quark model and extend our understanding of SU(3) symmetry breaking 
[ |l], 0. The discovery of the rare decay (ft — > 7/7 by CMD-2 [ || has brought 
the last piece into the otherwise complete mosaic of radiative dipole magnetic 
transitions between light vector and pseudoscalar mesons. This observation 
was later confirmed by the SND group [ Recently CMD-2 presented 
an improved measurement of the rate of the <fi — > 7/7 decay [ |J based on 
the total data sample accumulated in the 0-meson energy range. The first 
preliminary results by the KLOE collaboration on the Br(<f) — > 7/7) value 
reported recently [ || are consistent with those from CMD-2 and SND within 
the experimental uncertainties. 

Because of the rather low probability of the studied decay as well as 
serious background problems resulting in strict selection criteria, none of the 
above measurements have good statistical precision. One way to increase a 
sample of 7/7 events is to use other decay modes of the 7/ and rj mesons. In 
most of the previous papers [ 0, £|, ||, |6| events of the decay (p — > 7/7 were 
selected in the mode rj' — > 7T + 7T~rj, rj — > 77. The pure neutral decay chain 
rj' — > 7r 7r°?7, rj — > 77 has been also used by the KLOE group [ ||. In this 
paper we study the decay — > rj lr y with the CMD-2 detector at VEPP-2M 
in the channels: 

^ — > 7/7 , rj' — ► ?77r + 7r~ , rj — ► ix + ix~ix Q or (1) 
cf) — > 7/7 , rj' — ► ?77r + 7r~ , 77 — ► 7r + 7r~7. (2) 

The analysis is based on the data sample collected in three scans of the 
c.m.energy range from 0.984 to 1.060 GeV performed in winter 1997-1998. 
The total integrated luminosity of about 11.6 pb" 1 corresponds to approxi- 
mately 16 million meson decays. 

The general purpose detector CMD-2 has been described in detail else- 
where [0. It consists of a drift chamber (DC) [ |J and a proportional 



2 



Z-chamber [ |9|, both used for the trigger, and both inside a thin (0.4 X ) 
superconducting solenoid with a field of 1 T. 

The barrel calorimeter [ [Ttl which is placed outside the solenoid, consists 
of 892 Csl crystals of 6 x 6 x 15 cm 3 size and covers polar angles from 46° 
to 134°. The energy resolution for photons is about 9% in the energy range 
from 50 to 600 MeV. The angular resolution is about 0.02 radians. 

The end-cap calorimeter [ |ll| which is placed inside the solenoid, consists 
of 680 BGO crystals of 2.5 x 2.5 x 15 cm 3 size and covers forward-backward 
polar angles from 16° to 49° and from 131° to 164°. The energy and angular 
resolution varies from 8% to 4% and from 0.03 to 0.02 radians respectively 
for the photons in the energy range from 100 to 700 MeV. 

The luminosity was determined from Bhabha scattering events at large 
angles [ [12| . 



2 Data analysis 

The specific signature of the processes ([!]) and (Q) is the final state with four 
charged pions coming from the interaction region and two or more photons. 
It is well known that the meson decays do not provide such a signature 
with the exception of the following decays into kaons: 

(f) -> K° S K° L , K° s -> 7T+7T-, K\ -> tt+tt-tt and (3) 
-> K + K~, K* -> t^tt+tT, K T -> n T n°. (4) 

However, because of the large lifetime of kaons charged particles originating 
from the kaon decays have a broad distribution of the impact parameter 
with respect to the interaction point. It is this feature which will be used 
to suppress this background and control the contamination of the final event 
sample. 

Another background reaction is e + e~ — > ujtt , uo — > 7r + 7r _ 7r° followed 
by the Dalitz decay of one of the pions: 7r° — > e + e~ r ). In this case the 
distribution of the minimum space angle between the charged tracks has a 
characteristic peak at small angles. 

We selected events with four charged tracks coming from the interaction 
region: the maximum impact parameter of the tracks d = max.f =1 (r min .) is 
less than 1 cm and the vertex coordinate along the beam axis z vert is within 
±10 cm. For better reconstruction efficiency, tracks were also required to 
cross at least two superlayers of the drift chamber: |cos#| < 0.8. Two or 
more photons have to be detected in the calorimeter. For the selected events 
the kinematic fit was performed taking into account energy — momentum 
conservation and assuming all charged particles to be pions. 
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One of the main problems in the analysis is additional ( "fake" ) photons 
induced by the products of nuclear interactions of charged pions in the detec- 
tor material. The following method was used to suppress such fake photons. 
In the kinematic fit the energy resolution of photons was loosened to the 
value cte 1 = E y + 20 MeV. Thus, the photon energy was allowed to vary 
in a wide range during the fit. As the first step, the fit was performed in 
the hypothesis of the 7r + 7r~7r + 7r~77 final state. For events with more than 
two detected photons (N^ et > 2) the fit was repeated for all possible pairs 
of photons and the pair with the smallest xL-27 characterizing the fit quality 
was selected. Then, for these events, the assumption of the 7r + 7r~7r + 7r~777 
state was tested. For events with N^ et > 3, the combination of three photons 
with the smallest XL37 was se l ec ted. Finally, the signature 7r + 7r~7r + 7r~77 or 
7r + 7r~7r + 7r~777 was assigned to the event depending on whether the magni- 
tude of xL-27 or XL37 was sma U er - Then the requirement on the fit quality 
X 2 = min(x4 7r27 , xL-37) was a PPli e d ; X 2 / n -d-f < 10/4. Events with the re- 
constructed photon energy below 30 MeV were rejected from the subsequent 
analysis. 

The background coming from the decay — > K + K~ , in which products 
of kaon nuclear interaction scatter back to the drift chamber and induce 
two extra tracks or one of the kaon decays via the — > n ± ir + n~ channel, 
accompanied by fake photons, was suppressed using ionization losses of the 
tracks ( see [ |8| , [TBI for more detail). Requiring E^/2E}j earn > 0.7, where E ilT = 
^2i=i VPi + m2 r * s the total energy of the tracks, additional suppression of 
this type of the background was achieved. 

To suppress the background reaction e + e~ — > un , where uj — ► 7r + 7r~7r° 
and one of the neutral pions decays to e + e~7, a previously tested procedure 



for e/7r-separation (see [ 14, |Tf| for details) was employed. In this procedure, 
using the data samples of e + e _ — > 7r + 7r~7r° and e + e~ e + e~j events, the 
probability density functions for the parameter Ec s i/\p\ were obtained for 
both particle types (e/n) and signs (+/—): f w +, f n -, f e + and / e -. Here Ec s i 
is the energy deposition in the Csl calorimeter and p is the momentum of the 
particle with a track matching the cluster in Csl. Then the probability for the 
pair of tracks to be pions was defined: W n + % - = f n +fn-/(fn+fn- + fe+fe-)- 
Figure [l] shows the distribution of W^^- for the events of the processes 
e + e~ — ► e + e~7 and e + e~ — > 7r + 7r~7r° demonstrating clear separation of the 
two types of the processes. Since the decay ir° — > e + e~7 features a small space 
angle between e + e~, we searched for a pair of oppositely charged particles 
with a minimum space angle ip m i n . The probability W n + n - was calculated 
for this pair. After that, the condition W n + n - > 0.5 suppressed about 70% 
of the wit events. In addition, ionization losses dE/dx were used to reject 
con events with a low momentum of e + e~ [ [T3 . 



4 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 1: e/ir separation probability. Distribution of W / 7r + 7r - for the events of 
the processes e + e~ — > e + e~7 (histogram) and e + e~ — > tt + 7t~7t° (points with 
errors) 
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Figure 2: The invariant mass M^+^-. a) Points with error bars — data; 
transparent histogram — MC of the process (0); hatched histogram — MC of 
the reaction @; b) Simulation of the <p — > v'l decay; transparent histogram 
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About 250 events were selected after application of the criteria described 
above. For these events the distribution of the invariant mass of the pairs 
of charged pions M n + n - is shown in Fig.^a) by the black squares with error 
bars (four entries from one event). The transparent histogram in this Figure 
presents the distribution obtained by the complete Monte Carlo simulation 
(MC) of the CMD-2 detector [ |]| for the process (|), while the hatched 
histogram is the MC of the reaction (f|). The number of events in KgK® MC 
samples corresponds to approximately the same integrated luminosity as in 
our data. For comparison, in Fig.|b) the same distributions of M n + n - are 
shown for the simulation of the decays (Jj) and (Q). 

As one can see from Fig. ||]a), the main contribution to the selected events 
comes from the process <fi — > K S K® followed by the decays K s — > 7r + 7r~ and 
K\ — > 7T + 7r~7r°. The requirement that M n + n - < 450MeV/c 2 for each of the 
four possible tt + tt~ pairs efficiently rejects the K S K^ events and does not 
influence the events of the decay under study (see Fig. 0b)). After that only 
23 events survive. 
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Figure 3: The minimum space angle between the tracks versus the maximum 
impact parameter of the tracks. Black squares in both Figures present the 
experimental data; a) points are the MC of the decays (|l|) and (||); b) o is 
the MC of K° S K° L events; □ is the MC of K+K~ events; A is the MC of 
ujtt events. The number of events in the MC samples of charged and neutral 
kaons and the number of simulated uir events exceed the number of events 
of these reactions expected at our integrated luminosity by a factor of 10 and 
4 respectively 
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To complete separation of the 7/7 events and estimate the remaining 
background from the ujtt , KgK° and K + K~ events, the distributions of the 
minimum space angle between the tracks and the maximum impact parame- 
ter of the tracks were analyzed for 23 selected events. In Fig.^i) the scatter 
plot of ipmin versus d is shown for the experimental data in comparison to 
the simulation of the decays ([!]) and (|2|). In Figfjb) the same experimental 
data are compared to the MC of the background events KgK^ , K + K~ and 
ujtt . 

It is obvious from these Figures, that events of the decay under study 
contribute to the region of small impact parameters and feature a broad 
distribution of ipmin- The background events ujtt fall mainly in the region of 
small angles: ipmin ^ 0.2, while the events with neutral and charged kaons 
have approximately flat distribution of the maximum impact parameters of 
tracks with d > 0.3 cm. Correspondingly, the plane of the parameters ipmin 
and d was subdivided into three following regions: 

1. the 7/7 region: ip m in > 0.3 and d < 0.3 cm; 

2. the ujtt region: ipmin < 0.3 and d < 0.3 cm; 

3. the KK region: ipmin > 0.3 and d > 0.3 cm. 

The fourth region (ip m i n > 0.3 and d > 0.3 cm) contains a small contribution 
from all three background reactions. 

Ten candidate events of the <p — > rf decay were selected in the 7/7 region: 
N v ' y = 10. For these events the distribution of the missing mass for the 
pair of charged pions and the photon with the energy closest to that of the 
monochromatic photon from the decay <p — > 7/7 is shown in Fig. f|. The 
MC spectrum in this Figure was obtained as a sum of the contributions of 
the processes ([I]) and (0) with the numbers of events proportional to the 
corresponding decay probabilities. One can see that the experimental data 
agree well with the simulation. 

The number of background events was found to be: N^o = 3 and Nkk = 
8 in the regions ujtt and KK respectively. 

Figure |5|a) shows the distribution of the missing mass of all four charged 
pions versus the maximum impact parameter of the tracks for the events, 
which were selected under the condition ip m in > 0.3. Good agreement is ob- 
served between the Monte Carlo simulation of the signal and the experimental 
data for d < 0.3 cm. In Fig. |5|b the same distribution of the experimental 
data is shown in comparison to the simulation of the background processes 
(0) and (f|). These Figures provide one more verification of our assump- 
tion that the 7/7 events contribute to the region of small impact parameters, 
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while events with the decays of neutral and charged kaons have tracks which 
originate at rather large distances from the interaction point. 

The following probabilities to observe the background events ujti , K^K L 
and K + K~ in the 7/7 region, were obtained from the MC: 

w = iV M o(^>Q.3A0.3cm) = Q M ± Q 1Q 
N W K°KWmin < 0.3 , d < 0.3 cm) 

w N KR {%l) rnin > 0.3 , d < 0.3 cm) nin , nn , 
w KK — "77 — t, ; 7 = 0.10 ± 0.05 . 

N K K\Wmin > 0-3 , d > 0.3 cm) 

Using these probabilities and the numbers of events N wn o and N kR , the back- 
ground contribution into ten candidate events N^^ was estimated. Finally, 
the number of the 7/7 events was calculated from the formula: 

AV7 = ^Vr - N KK W KR - N^oW^o = 8.8±|:1 . 

The error is determined by the statistics of selected events only, as well as 
the background contribution whereas the errors of the probabilities W W7r o and 
Wifjf will be included in the systematic uncertainty. Here we did not take 
into account the possibility to observe 7/7 events in the background regions 
as well as the probabilities to detect the ujtt events in the region KK and 
vice versa. All these effects contribute less that 3% into N^^ and were also 
included into the systematic uncertainty. 

The detection efficiency was obtained using the Monte Carlo simulation 
of the processes (@) and (|): e 3w = 0.097 ± 0.003 and e^ 7 = 0.091 ± 0.003. 

The process @ was also used for the normalization and calculation of 
the — > 7/7 decay probability. To select KgK® events, conditions simi- 
lar to those for 7/7 were used with one exception. We required M m +^- > 
450MeV/c 2 for at least one 7T + 7r pair. Note that this condition effectively 
suppresses the process @. The selection efficiency was determined from MC: 
e K° s Kl — (4-9 ± 0.3) ■ 10~ 4 . The number of detected events is N K o K o = 216. 

The relative probability of the decay <fi — > 7/7 was calculated using the 
formula: 

Br((p -> 7/7) AV 7 Br(K° s -> ir + ir-)Br(K° L -> tt+tt-tt ) 



5r(0 K S K° L ) N K o K o BrU ^ tt+tt-t/) 

x !M ( 5) 

e 3lT Br(ri — > 7r+7r 7r°) + e nnJ Br(r] — > 7r+7r 7) 
= (1.46±g;|t ± 0.18) -10" 4 . 

The systematic uncertainty arises from the following contributions: 6.7% 
comes from the uncertainties in determination of the detection efficiencies e^, 
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Figure 4: The missing mass for the pair of pions and the photon with the 
energy closest to that of the monochromatic photon from the decay — > 7/7; 
histogram — Monte Carlo, black squares with error bars — data 
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Figure 5: The missing mass of Air versus the maximum impact parameter of 
the tracks. The black squares in both pictures - experimental data; (a) dots 
- sum of Monte Carlo of the processes (P and (b) light squares - MC of 
the <j) — > K Q S K° L decay, light triangles - MC of the <p — y K + K~ decay 
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Table 1: Summary of the (f) — > 7/7 decay rate measurements 



Experiment 


7/7 decay 
mode 


Br (<j>—tr)' y) 
Br(<f>— >r)7) ' 

10~ 3 


73 r (4>~*i7 '7) 

10" 4 


5r(0 — > 7/7), 
10~ 5 


CMD-2 [ | 


7/ — > 7T + 7T^77, 
77 -> 77 


9.5±^ ±1.4 




12.0^ ± 1.8 


SND [1 


77' — > 7T + 7T — ?7, 
7/ -> 77 






6.7^9 ±1.0 


CMD-2 [§, 
includes [ |3| 


7/' — > 7T + 7T^77, 
77 -> 77 


6.5i^±0.8 




8.2t££ ±1.1 


KLOE [§ 


7/' — > 7T + 7T — 77, 
77 -> 77 


7.1 ± 1.6 ±0.3 




8.9 ±2.0 ±0.6 


7/' — > 7T U 7T U 77, 
7/ -> 77 


6.81^ ± 0.9 






CMD-2, 
this work 


7/' — > 7T + 7T _ 77, 

77 — > 7r + 7r _ 7r or 

77 — > 7T + 7T _ 7 




1.46+^ ± 0.18 


4.9^ ±0.6 



e 7r7r7 and 7% comes from the stability of the selection criteria, 6.4% 

is due to the background subtraction and 4.1% comes from the uncertainties 
in the decay probabilities used in the Eq.(|5]). All decay probabilities were 
taken from [ [T7|]. The overall systematic uncertainty is 12.3%. 

Finally, taking the value of Br(<p —>■ K^K^) from [ 17], the following 
branching ratio was obtained: 

Br((f> -> rfi) = (4.9+^2 ± 0.6) • 10~ 5 . 

The energy behaviour of the cross section calculated using ten selected 
events does not contradict the shape expected for the sum of the meson 
and non-resonant background. 



3 Discussion 

The currently available data on the measurements of the <fi — > 7/7 decay 
probability are shown in Table [I] in comparison to our result. The value of 
the decay probability Br {(f) — > 7/7) obtained in our work is lower than that in 
but is consistent with their results within a statistical uncertainty. 
The previous CMD-2 measurement of Br [(f) — > 7/7) in the decay chain 
~ * v'li v' ~> n + -K~ri and 77 — > 77 [ [J is based on the data set, which is 
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statistically independent of our data sample. Analysis of the sources of sys- 
tematic errors in both measurements also shows that they are uncorrelated. 
Thus, to improve the accuracy of the branching ratio determination, we can 
combine both results: 

Br (<f) - r/ 7 ) C MD-2 = ( 6 - 4 ± L6 ) ' 10 ~ 5 ■ 

The authors are grateful to the staff of VEPP-2M for excellent perfor- 
mance of the collider, to all engineers and technicians who participated in 
the design, commissioning and operation of CMD-2. 
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